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Abstract 

Recently, the Columbia group obtained the valence quark mass dependence of the 
chiral condensate for lattice QCD simulations with dynamical fermions for a series of 
couplings close to the critical temperature. In this letter we show that the data for 
temperatures below Tc and sufficiently small valence quark masses can be rescaled to fall 
on a universal curve. Using chiral random matrix theory we obtain an exact analytical 
expression for the shape of this curve. We also discuss the universal scaling behavior of the 
mass dependence of the mass derivatives of the chiral condensate (scalar susceptibilities) . 
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Numerical simulations of QCD have contributed greatly to our present understanding 
of the QCD chiral phase transition (see for a recent review). It is now generally 
believed that the critical temperature is about 150 MeV, and that, for two light flavors, 
the phase transition is of second order ||^ . In the neighborhood of the critical temperature 
the particle spectrum changes drastically which may lead to important experimental 
implications. The hadron masses are closely related to the chiral condensate, the order 
parameter of the chiral phase transition, which, according to the Banks-Casher formula 

is proportional to the average spectral density of the Dirac operator in the vicinity of 



Unfortunately, the volumes of present day lattice calculations are relatively small. On 
the other hand, the smallest eigenvalues of the lattice QCD Dirac operator show the most 
significant volume dependence. In order to extract the chiral condensate reliably from 
lattice QCD simulations detailed knowledge of the volume dependence of the spectral 
density near zero is essential. Recent lattice calculations of the Columbia group 0] have 
probed this part of the spectrum via the so called valence quark mass dependence of the 
chiral condensate. In this paper we attempt to understand this dependence for masses 
well below the typical hadronic scale A. We also give results for the derivative of the 
chiral condensate (scalar susceptibilities) in this mass region. 

The valence quark mass dependence of the chiral condensate is defined by 



and is the total number of eigenvalues of the lattice QCD Dirac operator. The average 
over gauge field configurations, denoted by (■■■), is weighted by the gluonic action and 
the fermion determinant with sea quark mass equal to rrisea- It is clear that, for a finite 
number of eigenvalues, E^(m) ~ 1/m for asymptotically large valence quark masses. For 
valence quark masses well below the average position of the smallest eigenvalue, Amin, we 
expect that S^(m) ~ m. For masses larger than Amin but still sufficiently small we hope 
to see a plateau which can be identified as the chiral condensate (this assumes a smooth 



zero. 




(1) 



where p(A) is the average eigenvalue density of the Dirac operator 



p(A) = 5:(<5(A-A„)) 



(2) 
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sea quark mass dependence of the chiral condensate). The usual relation between the 
chiral condensate S and the average spectral density is given by the Banks-Casher 
formula 



(3) 



Here, p(0) is defined as the extrapolation to zero of the eigenvalue density many spacings 
away from zero. This can be done unambiguously if there are many eigenvalues below a 
typical hadronic scale A (a more rigorous definition is obtained by taking the limits m ^ 
and — s> oo in this order). Above this scale the eigenvalue density is dominated by the 
perturbative eigenvalue density ~ A^A^. In lattice calculations the latter dependence has 
to be adjusted for finite size effects. Because of the existence of a finite chiral condensate 
the Banks-Casher formula tells us that the smallest eigenvalues of the Dirac operator 
are spaced as ~ tt/A^S. This allows us to define a different limiting procedure of the 
eigenvalue density of the Dirac operator, namely, the microscopic spectral density defined 
by i 

1 z 

Ps(z) = lim p( ). (4) 

In terms of this density, eq. (1) can be rewritten as 

roo2miVSp5(^ 
^ ^""^ " Jo + m2Ar2E2 ' 

which suggests to plot the lattice data as S(m)/S versus mNH. 

In Fig. 1 the points show the original Columbia data obtained from a 16^ x 4 lattice 
with dynamical (staggered) fermions {Nf = 2) with mass mgeaO' = 0.01 (a is the lattice 
spacing). Results are given for three different values of /3 (see label of the figure) below 
Tc. The smallest eigenvalue is of order 2.5 10~^ for (3 = 5.245 and increases gradually for 
larger values of (3. A clear plateau is observed for two smallest values of (3. It is located 
around ma ~ 10"'^. The dotted curves are fits to the lattice data for masses larger than 
ma ~ 0.001. However, we did not fit directly the valence quark mass dependence but 
used a reasonable fitting function for the eigenvalue density, which allowed us to extract 

p(o)/iv. 
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Figure 1: The valence quark mass dependence of the chiral condensate. The points 
represent the Columbia data for values of /? as indicated in the label of the figure. Fits 
to their data are shown by the dotted curves. The result for the noninteracting Dirac 
operator is given by the dashed curve. 



In Fig. 2 we show the rescaled data for the same values of /3. Remarkably, they fall 
on a single curve suggesting that the microscopic spectral density is a universal function. 
Below we obtain an analytical expression for the shape of this curve (shown by the full 
line for N ^ = in this figure). 

From the analysis of spectra of complex systems we know that the correlations between 
eigenvalues of the Hamiltonian on a scale of no more than a finite number of level spacings 
can be described by a random matrix theory with the appropriate anti-unitary (time 
reversal) symmetry 0, §|, 0- The situation in the case of the Dirac operator is more 
complicated. First of all, because of the ?7A(l)-symmetry, in a chiral basis the massless 
Dirac operator has the block structure (for zero quark masses) 



T 
Tt 



(6) 

Since both the gauge fields and the Dirac matrices are in general complex, the anti-unitary 




Figure 2: The valence quark mass dependence of the chiral condensate E^(m) plotted as 
E^(m)/S versus mN'E. The dots and squares represent lattice results obtained by the 
Columbia group [1] for values of (3 as indicated in the label of the figure. The full lines 
correspond to the random matrix result (eq. (8)) for z/ = and Nf as indicated in the 
label of the figure. Also shown are results (dashed curves) for the quark mass dependence 
with equal valence and sea quark masses. 

symmetries act both in color and Dirac space. For a detailed discussion of this classifica- 



tion we refer to [|lOl- At this moment we only want to remark that for three colors with 
fundamental fermions there are no anti-unitary symmetries, leading to complex matrix 
elements Tjj. This ensemble was called the chiral gaussian unitary ensemble (chGUE). 



Because of the chiral structure of the Dirac operator, all nonzero eigenvalues occur in 
pairs ±A„. Therefore, A = is a special point in the spectrum. This allows us consider two 
types of level correlations: i) microscopic correlations in the bulk of the spectrum; they 
are given by the Dyson-Mehta-Wigner random matrix ensembles and ii) the microscopic 
spectral density (defined in (4)). At this moment it is not clear in how far deformations of 
the random matrix ensemble while keeping its chiral structure (6) lead to the destruction 
of both types of universal correlations at the same time, or that one of them is more 
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strongly universal than the other. Up to know we have not been able to construct a 



random matrix ensemble for which the latter is true |]13|. This is an important point 
because recently Kalkreuter obtained all eigenvalues of the Dirac operator for a 12^ 
lattice with dynamical (staggered) fermions. We have analyzed the eigenvalue correlations 
in the bulk of the spectrum and found [0 that they are given by the Dyson-Mehta-Wigner 
random matrix ensembles. If the two types of universal behavior are mutually inclusive, 
which we believe to be case, the microscopic spectral density of the Dirac operator is 
given by random matrix theory as well. Further arguments in favor of this conjecture 
can be found in [jl4[. General universality arguments show that the properties of 



random matrix ensembles are not sensitive to specific shape of the distribution of the 
matrix elements. For convenience, and in agreement with the maximum entropy principle 



Tql , we have chosen a gaussian distribution. 

The microscopic spectral density of the corresponding chiral Gaussian Unitary En- 
semble was obtained in [|l^, and is given by 

Ps{z) = -{Jn^+A^) - JNf+u+liz)JNf+u-liz)). (7) 

Here, Nj is the number of (massless) flavors and u is the topological charge. Using (5) 
we are able to calculate the valence quark mass dependence of the chiral condensate. The 
integrals are known analytically fl^ which leads to a remarkably simple result 



= u{lNj+u{u)KNf+u{u) + lNf+u+l{u)KMj+u-l{u)), (8) 

where the and the are modified Bessel functions and 

u = mNE. (9) 

The valence quark mass dependence for zero, one and two massless flavors and topological 
charge z/ = as given by eq. (8) is represented by the full lines in Fig. 2. Obviously, the 
lattice data are described by the curve with Nj = u = 0. In hindsight this is no surprise. 
First, in the region where the valence quark mass is much less than the sea quark mass, 
the fermion determinant has no bearing on the Dirac spectrum and we have effectively 
Nf = 0. Second, on a lattice the zero mode states and the much larger number of nonzero 
modes states are mixed, so that the effective topological charge that reproduces the lattice 



results is equal to zero. The small u behavior of the valence quark mass dependence is 
—u logM for Nf + = and linear for Nf + 1/^0. For valence quark masses smaller than 
the smallest eigenvalue of the Dirac operator over the finite ensemble of QCD gauge field 
configurations, a linear mass dependence is expected, which is indeed observed in lattice 
gauge calculations. 

For comparison we have also shown the sea quark mass dependence of the chiral 
condensate (dashed line in Fig. 2). This can be obtained most conveniently from the 



finite volume partition function which follows from general symmetry arguments [|18 
or can be derived starting from the chiral random matrix ensemble [|I9|. The effective 
partition function describes the mass dependence of the QCD partition function for space 
time volumes 

1 < yV4 ^ J_ (10) 
A m^r 

(A is a typical hadronic scale and ~ mA is the pion mass). It is straightforward to 
obtain the sea quark mass dependence of the chiral condensate from the effective partition 
function. In the sector of zero topological charge one finds |jl8 

fo. A^.^l. ,11) 

S Io{u) 

and 

These results are shown by the dashed curves in Fig. 2 and should be compared to lattice 
data for the sea quark mass dependence of the chiral condensate which, unfortunately, are 
not available at this moment. For m <^ Amin the mass in the determinant in not relevant 
and S'^('u) should approach S^('u) which is indeed seen in Fig. 2. 

For what range of masses can we expect the valence quark mass dependence to be 
described by random matrix theory? For Nf ^ 0, the small mass limit chiral random 
matrix theory reduces to the partition function of Leutwyler and Smilga ||19| , which gives 
the m-dependence of the QCD partition function in the range (10). The Nf = limit can 
be obtained from results at arbitrary Nf by taking the limit Nf — 0, in analogy with the 

replica trick in condensed matter physicsQ. Since the condition (10) is Nf independent 

^Alternative, one may use a supersymmetric formulation to cancel the fermion determinant leading 
to supersymmetric quenched chiral perturbation theory ||22|. 



it should be valid for Nf — as well (with m,r = "^A and m equal to the valence quark 
mass) . On the other hand, for lattice calculations to be physical, the pion should fit well 
inside the box, 

^ < L, (13) 



which leads to the condition 

m < < msea- (14) 

Because the pion is the hghtest hadron the condition yA <^ L in (10) is satisfied automat- 
ically. In practice the pion compton wave length is not much smaller than the length of 
the box. Indeed, the lattice data show a universal behavior up the valence quark masses 
close to the sea quark mass. 

It is natural to expect that also the derivatives of the chiral condensate with respect 
to the mass (the scalar susceptibilities) show a universal scaling. Below we will discuss 
the following susceptibihties 

1 1^11 1^1 

N ^ ' N jfj^^ iXk + miXi + m N ^^Xf. + 

1 1^1 

where the quark masses have been put equal to m after differentiation. The latter quantity 
is the isovector scalar susceptibility and x^^^ enters in the isoscalar scalar susceptibility, 
^ 2x^^" + x"""" (for two flavors). 

As follows immediately from (16) x^°^ can be calculated from average spectral den- 
sity (2). On the other hand, x^^^ can be expressed into the connected two-point level 
correlation function 

Pe(A,y) = (p(A)p(AO)-(p(A))(p(A')). (17) 
If the eigenvalues are uncorrelated we have 

p,(A,A') = (p(A))5(A-A'), (18) 

which leads to 

X'^^ = -. (19) 
m 
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In random matrix theory, the eigenvalues are strongly correlated, which gives rise to a 
quite different prediction. We will consider the susceptibilities in the limit mNJ] ^ 1, but 
with m <^ A. If, at the same time, m?Nh? <^ 1, the partition function of Leutwyler and 
Smilga or chiral random matrix theory can be used to calculate the susceptibilities ||20 



In this limit the leading order contribution is expected to be of the form 1/m^N. In the 
physical limit, where the m?NK^ correction to the vacuum energy cannot be neglected, 
general symmetry arguments are insufficient to determine the susceptibility. 

Since we are interested in the valence quark mass of the susceptibility we cannot use 
the effective partition function but have to rely on chiral random matrix theory. The 
valence quark mass dependence of x^°^ can be obtained from TX {u) given in eq. (8), 

X~»=S.„E>'(«)~^^i (20) 

where the asymptotic result is valid for u ^ 1 (the coefficient of the leading order term 
in lu vanishes). For Nf = 0, x'^'^ can be obtained from the two level correlation function 
| T5[| . The result is given by 



where both F'^'^(m) can be expressed in modified Bessel functions and approaches 1 for 
M > 1. 

In the limit m^NA^ ^ 1 the susceptibility is no longer determined by the lowest order 
terms in a chiral expansion of the vacuum energy. However, it is possible to calculate the 
susceptibilities in the framework of chiral random matrix theory but the results are model 
dependent, x™"^ can be obtained directly from the semicircular eigenvalue distribution. 
X"^^^ could only be obtained in an indirect way from the bosonized random matrix partition 
function. The results for all Nf are given by (for im?NK ^ 1) 

X^°° = -y, (22) 
x''^ = 0. (23) 

The result for x'^'^ is a consequence of the stiffness of the random matrix spectra. 

Unfortunately, no lattice data are available for the valence quark mass dependence of 
the susceptibilities quoted above. What has been calculated ||2^ is the dynamical quark 



mass dependence of the susceptibilities. Because, in this way only relatively large quark 
masses have been studied, ultraviolet divergent terms make direct comparison difficult. 
In rri^V ~ 1, but the numerical results show that, below the critical temperature. 



the quark mass dependence of the susceptibilities is ~ 1/m. One explanation is that the 
eigenvalues of the Dirac operator are only weakly correlated so that eq. (20) is applicable. 
On the other hand, spectral averaging shows strong correlations between the eigenvalues 
implying that the ensemble and the spectral average of the eigenvalue correlations 



are not equal. This exciting possibility deserves further attention. 

In conclusion, we have shown that the valence quark mass dependence of the chiral 
condensate is lattice QCD is given by a universal curve that can be described by chi- 
ral random matrix theory. As soon as lattice QCD results for the valence quark mass 
susceptibilities are available, we hope to extend our analysis to these observables as well. 
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